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LONG-TERM  GOALS 

The  general  goal  of  this  project  is  to  examine  the  seabed,  quantitatively  document  the  fluxes  of  fine 
sediment  (over  different  time  scales),  and  thereby  validate  localized  measurements  and  numerical 
models  of  sediment  transport. 

OBJECTIVES 

The  specific  objectives  are  to: 

a)  document  changes  in  bed  elevation  (deposition,  erosion)  on  time  scales  intrinsic  to  the  driving 
forces;  e.g.,  tidal  cyclicities,  wind-driven  currents  and  local  waves,  river  discharge,  and 
interannual  variability; 

b)  measure  net  accumulation  rates  over  decades  at  many  sites  to  calculate  fine-sediment  budgets  for 
both  Willapa  and  Skagit  tidal  systems; 

c)  examine  sedimentation  at  sufficient  locations  to  characterize  spatial  variability  of  grain  size  and  its 
vertical  stratification. 

APPROACH 

The  tidal-flat  sedimentation  in  Willapa  mud  flats  and  Skagit  sand  flats  will  be  contrasted,  with  a  focus 
on  understanding  the  import  of  mud  to  the  Willapa  flats  and  the  export  of  mud  from  the  Skagit  flats. 
This  involves  development  of  sediment  budgets  to  document  quantitatively  the  fate  of  muddy 
sediment,  and  also  includes  investigation  about  the  variability  for  sedimentation  of  several  time  scales. 

Because  UW  is  located  near  the  Willapa  and  Skagit  tidal  flats,  we  are  helping  to  provide  logistical 
support  to  colleagues  traveling  long  distances  to  work  in  these  areas.  We  have  also  helped  obtain  the 
substantial  number  of  permissions  needed  for  the  research,  and  have  acted  as  the  primary  contact  for 
the  federal,  state  and  county  regulatory  agencies. 
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WORK  COMPLETED 


During  2008-2009,  three  extended  surveys  were  completed  in  Willapa  Bay  for  detailed  observations 
and  sampling:  Sep  2008,  one  week  (spring  tide);  March  2009,  two  weeks  (spring  and  neap  tides);  July 
2009,  two  weeks  (spring  and  neap  tides).  In  addition,  regular  trips  were  taken  to  Willapa  Bay 
(approximately  once  per  month)  for  time-series  sampling  and  for  testing  newly  constructed  equipment. 
As  part  of  the  project,  a  mobile  barge  was  designed  and  built  for  spatial  sampling  in  shallow  water,  and 
also  for  semi-diurnal  time-series  observations  at  fixed  locations  (using  its  jack-up  capability. 
Additionally,  a  coring  device  was  constructed  that  allowed  collection  of  cores  (40  cm  long)  from  boat 
or  barge,  with  little  impact  on  the  muddy  tidal  flat.  The  cores  are  easily  examined  by  x-radiography 
and  subsampled  for  radiochemical  and  grain-size  analyses. 

A  number  of  factors  necessitate  that  the  sampling  strategy  for  the  Skagit  Bay  tidal-flat  system  be 
different  than  in  Willapa  Bay.  A  major  distinction  is  the  direct  input  from  the  Skagit  River  of  much 
new  sediment  (3-4  million  tons/year),  which  is  supplied  by  a  bimodal  hydrograph.  During  the  late 
spring  and  early  summer,  gradual  snow  melt  in  high  elevations  provides  much  freshwater  and  fluvial 
sediment.  During  the  late  autumn  and  early  winter,  rainfall  events  (at  times  coupled  with  rapid  snow 
melt  in  low  elevations)  provide  episodic  floods  of  the  river  with  much  suspended  sediment.  Small 
boats  were  used  to  sample  the  tidal  flat,  and  sampling  occurred  during  these  two  periods  of  high 
discharge:  Jun-Jul  2008,  Nov  2008-Jan  2009,  and  Jun-Jul  2009.  In  addition,  cruises  on  ships  were 
undertaken  to  the  subtidal  areas  of  the  dispersal  system  (e.g.,  Saratoga  Passage,  Deception  Pass)  in 
order  to  document  Skagit  sediment  that  bypassed  the  tidal  flats.  These  cruises  occurred  during  the 
periods  between  high  discharge:  Aug  2008  (R/V  Barnes),  Feb  2009  (R/V  Centennial),  Aug  2009  (R/V 
Centennial)  and  on  two  academic  cruises  in  Apr  2008  and  Apr  2009  (both  on  the  R/V  Bames). 

RESULTS 

Willapa  Bay 

Observations  have  documented  tidal  asymmetry  of  sediment  transport  on  the  Willapa  flats.  For 
example,  during  the  latter  stages  of  ebbing  tide  in  March,  much  mud  was  moved  into  and  down 
channels  as  high-concentration  flows.  The  concentrations  of  suspended  sediment  were  much  reduced 
in  July.  During  March,  x-radiographs  in  channel  thalwegs  reveal  distinct  physical  stratification  with 
laminae  1-3  cm  thick,  which  were  likely  fonned  during  the  ebb  flows. 

On  seasonal  time  scales,  evidence  of  variability  was  observed  in  several  ways.  The  floors  of  secondary 
tidal  channels  (e.g.,  tributaries  to  Bear  Channel)  showed  a  dramatic  change  between  winter  conditions 
(March,  2009)  and  summer  conditions  (July,  2009).  The  muddy  channel  bottoms  with  physical 
stratification  (laminae  1-3  cm  thick)  gave  way  to  a  shell-hash  pavement  (Fig.  1)  in  summer.  The 
presence  of  a  shell-hash  layer  ~20  cm  below  the  winter  channel  bed  suggests  that  seasonal  input  and 
removal  of  sediment  (tens  of  centimeters  thick)  is  a  regular  occurrence. 
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Fig.  1  -  C  channel  of  Willapa  Bay  during  July  2009,  with  shell  hash  at  base. 

This  can  be  seen  in  top  of  x-radiograph  shown  to  left,  and  is  similar  to  the  buried  shell-hash 

layer  observed  during  winter. 

7Be  (half  life  53  days)  is  an  indicator  of  recent  sediment  input  from  a  terrestrial  source.  7Be  was 
commonly  found  in  the  surficial  seabed  after  winter  flooding  (i.e.,  samples  collected  in  March),  and 
indicated  import  of  sediment  likely  from  local  streams.  Th  (half  life  24  days)  on  particles  is  an 
indicator  of  scavenging  from  ocean  water  (with  dissolved  234Th)  and/or  sediment  that  had  been  in  the 
ocean  outside  Willapa  Bay  soon  before  collection  of  seabed  samples  (<3  mos).  Th  was  found  in 
almost  all  surficial  sediment  samples  collected  during  winter  2009. 

Numerous  cores  (~40)  have  been  collected  for  measurement  of  sediment  accumulation  rates  by  210Pb 
geochronology  (half  life  22.3  y).  Eventually  these  will  be  used  to  create  a  budget  for  sediment 
accumulating  on  the  tidal  flat  in  southern  Willapa  Bay.  The  cores  have  been  collected  with  clear 
recognition  of  their  location  with  respect  to  local  morphology  (e.g.,  channel,  bank,  flat),  due  to 
expected  differences  in  sedimentary  processes. 

Several  signatures  have  been  recognized  so  far.  These  include  cores  that  seem  to  have  broken  records 
of  accumulation  -  steady-state  sediment  accumulation  for  20-40  years,  following  an  erosional  break 
(Fig.  2).  Accumulation  rates  are  typically  3-7  mm/y,  although  rates  <2  mm/y  are  found  on  some 
portions  of  the  flats  and  some  sites  in  channels  (e.g.,  cut  banks)  are  erosional.  Accumulation  rates  in 
other  areas  of  the  flats  seem  to  be  very  rapid.  Nearly  uniform  excess  210Pb  is  observed  to  extend  deep 
into  the  seabed  (>40  cm),  which  makes  accumulation  rates  difficult  to  calculate  because  present  cores 
are  too  short.  But  the  thickness  of  excess  210Pb  likely  indicates  high  rates  of  sediment  accumulation 
(>2  cm/y)  in  these  places.  Vibracores  have  been  collected  and  will  provide  longer  210Pb  profiles. 
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Fig.  2  -  Core  from  Willapa  central  channel.  Indicates  accumulation  at  a  rate  of  4.5  mm/y  for  past 
~40y.  Break  in  Pb  profile  at  ~20  cm  depth  may  indicate  this  accumulation  followed  a  period 
of  erosion.  Dominance  of physical  stratification  in  x-radiograph  precludes  significant 

impact  of  bioturbation. 

All  cores  collected  for  measurement  of  sediment  accumulation  rates  have  been  examined  by  x- 
radiography,  and  typical  stratification  is  demonstrated  by  Fig.  2.  Physical  stratification  dominates  the 
channels,  banks  and  most  of  the  flat  surface.  However,  some  parts  of  the  flats  show  impact  from 
bioturbation,  and  are  mottled  or  homogenous.  The  sediment  grain  size  is  predominantly  silt  and  clay, 
but  the  distinct  physical  laminae  indicate  that  the  grain  size  of  deposited  sediment  changes  through 
time.  Despite  the  dominance  of  silt  and  clay,  samples  from  the  north-central  portion  of  the  flat  contain 
some  sand.  Future  grain-size  analyses  will  quantify  sediment  texture  at  all  core  sites. 

Skagit  Bay 

The  surface  of  the  Skagit  tidal  flats  commonly  has  ripples  associated  with  tidal  currents  and/or  surface 
waves.  These  features  demonstrate  the  dynamic  and  sandy  character  of  the  flats.  Following  high 
discharge  events,  ephemeral  mud  deposits  were  found  several  centimeters  thick  and  containing  7Be.  In 
addition,  graded  sand  layers  (i.e.,  upward  fining)  were  found  ~10  cm  thick,  suggesting  input  from 
suspension.  The  ephemeral  mud  and  graded  sand  deposits  were  located  along  the  flanks  of  the  Skagit 
North  Fork. 

7Be  is  typically  associated  with  silt  and  clay  particles,  because  they  are  more  chemically  reactive  than 
sand  particles.  The  distribution  of  7Be  following  discharge  events  shows  that  most  silt  and  clay 
bypasses  the  sandy  portion  of  the  Skagit  flats,  and  is  especially  common  on  the  outer  fringe  of  the  flats. 
Surface  samples  collected  on  the  winter  cruise  following  high  discharge  events  also  revealed  7Be  in 
northern  Saratoga  Passage  and  west  of  Deception  Pass. 
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Despite  the  ephemeral  muds,  no  thick  or  extensive  mud  deposits  have  been  observed  on  the  tidal  flats 
at  the  river  mouth.  Muds  are  transported  to  the  subtidal  region,  and  accumulate  there.  The  only 
subtidal  portion  of  Skagit  Bay  is  a  narrow  channel  along  the  west  side  of  the  bay,  near  Whidbey  Island, 
but  extensive  deposits  of  modern  mud  are  not  found  there  either.  Much  mud  is  transported  into 
Saratoga  Passage,  and  is  accumulating  there  at  rates  ~1  cm/y  in  the  northern  portion  of  the  passage 
(Fig.  3).  Excess  Pb  is  also  found  in  muddy  patches  west  of  Deception  Pass,  suggesting  that  some 
Skagit  fine-grained  sediment  is  leaving  in  that  direction. 
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Fig.  3  -  Fasten  core  from  northern  Saratoga  Passage.  Total  Pb  activity  decreases  to  supported 
levels  at  ~100  cm  depth  in  core.  Excess  activity  indicates  an  accumulation  rate  of  ~1  cm/y.  The 
seabed  is  ~60%  clay  and  40%  silt,  and  benthos  are  able  to  bioturbate  intensely  and  destroy  physical 
sedimentary  structure  (as  seen  in  X-radiograph  negative). 

The  flats  are  primarily  composed  of  sand,  with  fine  and  very  fine  sands  dominating  (2-4  phi,  64-256 
pm;  Fig.  4).  The  sediment  buried  in  the  flat  is  dominated  by  cross  bedding,  with  very  little  evidence  of 
bioturbation.  Mud  within  the  seabed  is  largely  restricted  to  thin  layers  (few  centimeters  thick),  which 
is  surprising  because  the  Skagit  River  discharges  millions  of  tons  of  mud  each  year.  Most  of  the 
sediment  in  Saratoga  Passage  is  mud  (Fig.  3)  and  this  area  is  the  likely  destination  for  much  of  the 
missing  sediment,  possibly  along  with  the  region  outside  Deception  Pass. 

IMPACT/APPLICATIONS 

The  research  completed  in  this  project  leads  to  an  improved  understanding  of  the  processes  that  control 
the  geometry  of  sedimentary  deposits  over  multiple  time  scales  and  in  diverse  oceanographic  settings. 
The  tidal-flat  studies  will  provide  an  understanding  of  rate  and  type  of  sediment  (mud,  sand)  emplaced 
on  tidal  flats,  which  are  critical  for  understanding  the  surface  morphology  (e.g.,  channel  configuration) 
and  internal  character  (e.g.,  sediment  stratification  that  affects  strength). 
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Fig.  4  -  Grain-size  distributions  for  two  typical  surface  samples  from  Skagit  tidal  flats. 
Lower  sample  (953)  is  well  sorted  fine  sand.  Upper  sample  (946)  has  a  tail  extending  into 

very  fine  sand  and  silt. 


TRANSITIONS 

Other  investigators  in  the  Tidal  Flats  DRI  are  transferring  the  results  from  this  effort  into  their  projects. 
Those  studying  the  seabed  incorporate  radiochemical  and  textural  data  to  document  the  processes  (e.g., 
physical  reworking,  bioturbation)  impacting  seabed  characteristics.  Researchers  analyzing  boundary- 
layer  processes  also  utilize  these  data  to  determine  instrumentation  sites.  Accumulation  rates, 
sediment  budgets,  and  grain-size  data  are  key  components  to  input  parameters  for  numerical  models. 

RELATED  PROJECTS 

Related  projects  include  studies  of:  the  seabed  by  R.  Wheatcroft  and  P.  Wiberg;  boundary-layer 
processes  by  A.  Ogston,  R.  Geyer,  P.  Traykovski,  and  D.  Ralston;  suspended-sediment  dynamics  by  P. 
Hill,  B.  Law  and  T.  Milligan;  seabed  thennal  processes  by  J.  Thomson  and  C.  Chickadel. 
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